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1 Introduction 

Demand for the multiloop calculations is constantly growing nowadays. As the num- 
ber of loops increases, the calculations become more and more complicated and 
require, almost necessarily, some stages to be done automatically using computers. 
The reduction of the loop integrals with the help of the integration-by-parts (IBP) 
identities [1, 2] is an important stage which can be automatized. This reduction 
allows one to reduce the calculation of the multiloop diagrams to the calculation of 
a finite set of master integrals. It is important that the reduction also allows one to 
obtain differential and difference equations, which can be used for the calculation of 
the loop integrals without explicit itegration. 

One of the most successful methods of the IBP reduction is the Laporta algorithm 
[3]. The algorithm is easy to implement and to use, given a sufficient amount of time, 
it works flawlessly. It also allows for a number of programming improvements. These 
advantages explain why many modern most powerful reduction programs heavily rely 
on this algorithm, in particular, AIR [4], FIRE [5], Reduze [6, 7], and many private 
versions. However, there are some weak points of this algorithm, which may put some 
restrictions on its application. First, the reduction generates heavy-weight databases 
of the discovered rules which can be too expensive to save. Therefore, typically, 
the reduction is performed each time "from the beginning" which requires the same 
identities to be solved in each run. Another disadvantage of the Laporta algorithm 
is due to the huge redundancy of the IBP identities [8]. This redundancy results in 
many unnecessary calculations which make the reduction slow. 

Another approach to the reduction is a derivation of the symbolic rules which 
can be applied to any problem of a given class. Its advantages are obvious: nothing 
is being solved in the process of reduction, therefore, the reduction is very fast. 
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Symbolic rules are small in size, so, they can be easily saved for future calculations. 
The bottleneck of the approach is the search of the symbolic rules, a stage which 
seems to require a lot of manual work. 

An ideal solution of the reduction problem would be, therefore, designing a 
program which automatically finds the symbolic reduction rules at the first stage 
and then uses those rules for the reduction. Much effort has been devoted to the 
developement of the approach connected with the notion of the Groebner basis [9- 
12]. This is due to the fact that the problem of IBP reduction is very similar to 
the reduction of the elements of some algebra with respect to its ideal. In the 
latter problem there is a known algorithm of the construction of the Groebner basis 
— the Buchberger's algorithm and its generalizations. When the Groebner basis 
is constructed, the reduction can be performed unambiguously and very fast. If 
this approach always worked for the IBP reduction, it would definitely be the most 
complete solution of the reduction problem. However, there is a small peculiarity 
of the IBP identities which results in the essential obstacle when generalizing the 
Buchberger's algorithm to the IBP reduction. As it was shown in Ref. [8], the IBP 
reduction problem can be reduced to the problem of reduction of some noncommuting 
polynomial ring with respect to the direct sum of the left and right ideals. While 
the Buchberger's algorithm seperately works for both left and right ideals, it appears 
to be difficult to generalize this algorithm to the desired case. Probably, the only, 
partly successful, attempt of this approach has been made in FIRE, where the notion 
of s-bases [11, 12] have been used. 

Of course, there are many peculiar features of the IBP identities which make the 
IBP reduction not the general-case reduction with respect to the direct sum of the left 
and right ideals. In particular, as it was shown in Ref. [8], the generating set of the 
left ideal appears to be equipped with a Lie-algebraic structure, and that of the right 
ideal consists of the commuting elements. It is quite possible that at some point in 
the future a general solution of the IBP reduction problem will appear. Meanwhile, 
one can try to develop some heuristic agorithms which are not guaranteed to work 
for each case, but, nevertheless, are useful from the practical point of view. With 
the lack of a systematic approach (i.e., a strict algorithm), this developement can be 
quite challenging. 

This short note describes a Mathematica package LiteRed which can be con- 
sidered as an attempt of the implementation of the heuristic approach to the IBP 
reduction. The package can be downloaded as a zip archive from 

http : //www. inp.nsk.su/~lee/programs/LiteRed/ 

2 General setup 

Assume that we are interested in the calculation of the L-loop integral depending on 
the E external momenta pi, . . . ,pe- There are N = L(L + l)/2 + LE scalar products 
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depending on the loop momenta If. 



h ■ Qj , 1 ^ i ^ L,j ^ L + E, 



(2.1) 




N 



(2.2) 



where the scalar functions are linear polynomials with respect to Sij. The func- 
tions Da are assumed to be linearly independent and to form a complete basis in the 
sense that any non-zero linear combination of them depends on the loop momenta, 
and any Sik can be expressed in terms of D^- Thus, each integral is associated with 
a point in Z^. Some of the functions correspond to the denominators of the 
propagators, the other correspond to the irreducible numerators. E.g., the i^-legged 
L-loop diagram corresponds to E = K — 1 and the maximal number of denominators 
is M = E + 3L - 2, so that the rest N - M = (L - 1)(L -f- 2E - 4)/2 functions 
correspond to irreducible numerators. 

IBP identities The IBP identities [1, 2] are based on the fact that, in the dimen- 
sional regularization, the integral of the total derivative is zero. They are derived 
from the identity 



Performing the differentiation on the right-hand side and expressing the scalar prod- 
ucts via Da, we obtain the recurrence relation for the function J. 

LI identities There is also another class of identities, called Lorentz-invariance 
(LI) identities due to the fact that the integral (2.2) is Lorentz scalar [13]. They 
have the form 



The differential operator in braces is nothing but the generator of the Lorentz 
transformation in the linear space of scalar functions depending on pk. Again, per- 
forming the differentiation on the right-hand side, we obtain LI identity. Though 
these identities can be represented as a linear combination of the IBP identities [8], 
they prove to be useful in real-life reduction. 

Sectors The notion of sectors can be introduced as follows. The 6 = {9i, . . . ,9]\r) 
sector, where 9i = 0, 1, is a set of all points (rii, . . . ,12^) in whose coordinates 
obey the condition 
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In particular, the point (6'i, . . . , Ojy) belongs to the {9i, . . . , 9^) sector, and will be 
referred to as the corner point of the sector. Owing to this definition, the integrals 
of the same sector have the same set of denominators. 

Scaleless integrals The scaleless integral can be defined as the one which gains 
additional non-unity factor under some linear transformation of the loop momenta. 
Obviously, if j (6'i, . . . , 9n) is scaleless, then all integrals of the sector (^i, . . . , 9n) are 
zero. We will call such a sector a zero sector. A simple and convenient criterion of 
zero sectors has been formulated in Ref. [8]. According to this criterion, the sector is 
zero if the solution of the IBP equations in the corner point (^i, . . . , 9]\f) result in the 
identity j {9i, . . . , 9^) = 0. Note that this criterion may miss some scaleless sectors. 
Though this seems to be a small problem (undetected zero sectors are simply reduced 
to lower sectors), let us explain on a simple example why this happens. 
Let us consider the massless one-loop onshell propagator integral 



Obviously, this integral is zero for any ui and n2. However, it can be explicitely 
checked that there is no linear combination of the operators di ■ I and di ■ k which 
acts as identity on the integrand for ni = ^2 = 1. In other words, the solution of 
the IBP identities in the corner point of the sector (1, 1) does not result directly to 
J(l,l) = (though, it results to, e.g., J(l,l) oc J (0, 2) ). In order to prove that 
the integral J (1, 1) is scaleless, let us consider instead the following operator 



where k is an auxiliary vector chosen to satisfy the conditions k"^ = and k-k = 1. It 
is easy to check that Oj (1, 1) = {d — 4) j (1, 1). Since the operator O is a generator 
of the linear transformation I I + e (^l + {I ■ k) k — (^l ■ k^ k^, the integral j (1, 1) is 
scaleless. The reason why the IBP identities failed to lead to the identity J (1, 1) = 
is that the construction of this identity required introduction of the auxiliary vector 
k. There is an interesting open question: whether the introduction of the auxiliary 
vectors (or tensors) may lead to a new kind of the identities independent of IBP 
identities (unlikely) or dependent on them, but still useful for real-life reduction. 

Symmetry relations In many cases there exist nontrivial linear transformations 
of the loop momenta which map the set of the denominators of one sector on the set 
of the denominators of the same, or another, sector. Those transformations have the 
form 



where Mij is a L x [L + E) matrix. It is easy to understand that, for nonzero 
sectors, there is only a finite number of such transformations, all subjected to the 
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condition 



= 1. Those transformations induce some mappings of 
the denominator set and also mappings of the numerators into the hnear polynomials 
of Di . These mappings give nontrivial identities between the integrals of two different 
(or one) sectors, which are conventionally called symmetry relations (SR). 

Ordering For the reduction procedure to work, it is necessary to define some suit- 
able ordering of the integrals, i.e., the ordering in Z^. It is natural to consider the 
integrals with less denominators to be simpler. This defines a partial ordering of the 
sectors. We can extend this ordering to the complete one by, e.g., saying, that the 
two sectors 6i and O2 with equal number of denominators are ordered lexicograph- 
ically. For the integrals in sector = (6*1, . . . , 9^) with K = '^^=1 ® ('^q ~" l/^) 
denominators we choose the folowing ordering. For the integral J (ni, . . . , njy) we 
determine the ordering weight w = Wq, Wi, . . . , w^), where w_i = J2a=i 

is a total power of the denominators and numerators, wq = X]q=i (1 ^ ^q) is the 
total power of the numerators, and {wi, . . . , wn) is obtained from (— |ni| , . . . , — IunD 
by shifting powers, corresponding to denominators, to K left-most positions. Then, 
the two integrals in one sector are compared by comparing the left-most dictinct 
entries of their ordering weights. In particular, the simplest integral in the sector 

(01,..., M is J(^^i,...,M- 

Of course, the choice of the ordering is not unique. In principle, for the possibility 

of the reduction, there is only one strictly required property of the ordering. It is 

that for any integral J (n) there is only finite number of the integrals simpler than 

it (this number depends, of course, on n). However, the following condition of the 

chosen ordering is essential for our consideration. For a given sector all components 

of the ordering weight w are just linear combinations of n^. In particular, it leads to 

the fact that the relation J (rii) -< J {112) between the integrals in the same sector 

is invariant with respect to the shift of both rii and 112 by some 6n, provided that 

J (ni + 6n) and J {112 + Sn) also belong to the same sector. 

Differential equations As it was mentioned above, the differential equations can 
be used for finding the master integrals. The simplest type of such equations is the 
differential equation with respect to the mass. Probably, the first example of their 
application is presented in Refs. [14-16]. The differential equations with respect to 
the invariant constructed of the external momenta have been introduced and applied 
in Refs. [13, 17, 18]. The peculiarity of the latter case is due to the fact that, though 
the integral depends on the external momenta only via their scalar products, the 
integrand also depends on the scalar products of the external momenta and loop 
momenta. Therefore, before differentiating the integral, it is necessary to express 
the derivative with respect to the invariant via the derivatives with respect to the 
external momenta. For example, if the integral depends on p"^, and p ■ q, the 
derivative with respect to p ■ g at fixed and q^ can be expressed in two equivalent 
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ways 



-J in) = — ■ — J n = — ■ — J n 



d{jp-q) {p.qf p2q2 dp [p . qf - p2q2 dq ' 

In general case, when there are E > 2 external vectors, we have the following formu- 
las: 

^ J{n) = J2 [G-']^^p.-d,,J{n) = J2 [G-']^^p.-d,,J{n) , 



d {pi ■ p2 
d 



J(n) = i5^[G-^]^^p,-9p,J(n) . (2.6) 



d{pl) ' ' 2 

/ p\ ■■■ pi- Pe\ 



where G = G (pi, . . . ,pe) 



is a Gram matrix. 



\pi-Pe ■■■ P% J 

Acting by the operator on the right-hand side on the integrand and performing 
the IBP reduction, one obtains the differential equation for t/(n). 

Dimensional recurrences Dimensional recurrences have been introduced in Ref. 
[19]. Since their introduction, they have been successfully applied to the calcultion 
of the different integrals, see Refs. [19-21]. Recently, a method of calculation of 
the multiloop master integrals, based on the dimensional recurrences and analytical 
properties of the multiloop integrals as functions of ci, has been introduced in Ref. 
[22] and successfully applied in Refs. [23-30]. For further purposes it is convenient 
to introduce the operators and B^, see Ref. [8], acting as follows 

(Aj/^)) (ni, . . . , riiv) = rij/^) (rii, . . . , rij + 1, . . . , ^at) , 

(5, J(^)) (ni, . . . , n^) = J(^) (m, . . . , - 1, . . . , n^) . (2.7) 

The original derivation of the dimensional recurrence relation in Ref. [19] relied 
on the parametric representation. The result of this derivation for the integrals which 
can be presented as a graph has a nice and compact form 



j('^-2)(n)=/.^^(A,...A,,jW)(n), (2.8) 



trees 



where ii, . . . ,ii numerate the chords of the tree, and fi = ±1 for the Euclidean/Minkovskian 
case, respectively. For computer implementations, probably, a more convenient for- 
mula has been derived in Ref. [26] . It reads 

J^'-'^ (n) = (/x/2)Met |2^-|^^|,,=i.„^| (n) . (2.9) 

Note that this formula is valid also for the integrals with numerators and dots. We 
will call Eqs. (2.8) and (2.9) the raising dimensional recurrence relations. 
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For completeness, we present also the lowering dimensional recurrence relation, 
obtained in Ref. [22] 

J''"-'' (n) = ^^%^^E-L+tr ^ ■ ■ ■ ' ^""^ ' ^'-'^^ 

where = a (a + 1) ... (a + -L — 1) is the Pochhammer symbol, V (f i, . . . Vk) = 
det G (f 1, . . . f fc) is the Gram determinant, and P {Di, . . . , Dj^) = V{qi, . . . qi+E)- 

If we consider one of these dimensional recurrence relations for the master integral 
and make the IBP reduction of the right-hand side, we will obtain the difference 
equation for this master integral. 



3 Using the LiteRed package 

A typical package usage includes two stages: the search of the reduction rules 
and their application for the reduction. During the first stage the definitions, related 
to the basis, in particular, the reduction rules, can be saved to the disk. These 
definitions should be loaded on the second stage and applied for the reduction. The 
basic example of a program searching the reduction rules is presented in Fig. 1. This 
program finds the reduction rules for the two-loop onshell massive propagator. 

Let us provide some comments. At each stage the program generates some ob- 
jects which are then used at the later stages. E.g., the command 
NewBasis [p2 , {sp [p-1] -1 , sp [p-l-r] -1 , sp [p-r] -1 , 1 , r} , {l , r}] sets up a new ba- 
sis p2, consisting of the functions Da depending linearly on the scalar products in- 
volving loop momenta 1, r (sp [a,b] stands for the scalar product of a and b, sp [a] 
is a shortcut for sp [a, a] ). This procedure checks that the set of is linearly inde- 
pendent and complete and generates several objects Ds [p2] , SPs [p2] , LMs [p2] , 
EMs [p2] , To j [p2] . The meaning of these objects should be clear from the output 
immediately following the NewBasis command. The objects IBP[p2] and LI [p2] 
generated by GeneratelBP call give the functions which return the Integration-By- 
Parts identities and Lorentz-Invariance identities in a given point, see Eqs. (2.3) 
and (2.4). E.g., IBP[p2] [nl,n2,n3,n4,n5] gives a list of IBP identities in a general 
point. The IBP identities generated by GeneratelBP procedure are necessary for the 
detemination of zero sectors by the AnalyzeSectors procedure. This procedure gen- 
erates the list ZeroSectors [p2] of zero sectors as well as the list SimpleSectors [p2] 
of the simplest nonzero sectors. The latter is used in the procedure FindSymmetries. 
This procedure finds internal and mutual symmetries of the nonzero sectors and gen- 
erates the lists MappedSectors [p2] and UniqueSectors [p2] . Any sector from the 
former list can be mapped onto some sector from the latter. The substitution rules 
for such a mapping can be found in j Rules [p2, . . .] for each sector js[p2, . . .] 
from the list MappedSectors [p2] . 
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The central procedure which tries to construct the complete set of symbolic 
rules for a given sector is Solve j Sector which we apply to each unique sector in 
this example. The result of its work is the set of rules for each js[p2, . . .] from 
the UniqueSectors [p2] list. All rules found are saved in the directory "p2 dir" and 
ready to use for the reduction. The typical program performing the reduction is 
shown in Fig. 2. 

Drawing graphs There is a possibility to draw graphs, corresponding to the 
integrals and sectors. We remark that, in principle, the graph is determined, up 
to some equivalences, by the set of internal lines. The possibility to automatically 
determine the graph, corresponding to the set of denominators, is planned in the 
future versions of the package. Meanwhile, the present version implements the fol- 
lowing. After defining the basis, one can attach a graph to the highest sector(s) 
(which can be depicted as a graph) by the command AttachGraph. Then the graph 
for all subsectors is determined automatically. 

For the above example there are no irreducible numerators, so the sector 
j s [p2 ,1,1,1,1,1] can be depicted by the graph. Then, the graph is attached by 
the command 

AttachGraph[js[p2, 1, 1, 1, 1,1], { 
(*5 internal lines*) 
{1^2, "I"}, 
{2^ 3, "1"), 
{3^ 4, "I"}, 
{1 ^ 3, "0"}, 
{2 -> 4, "0"}, 
(*2 external lines*) 
{0^1, "p"), 
{4^0, "p"} 

}]; 

Then, a graph of, say, sector j s [p2 ,1,1,1,1,0], can be drawn with the com- 
mand GraphPlot [ j Graph [ j s [p2 ,1,1,1,1,0]]]. Note that GraphPlot is a standard 
Mathematica function and the presentation of the graph can be altered using its op- 
tions, as can be found in the examples distributed with the package. 

Additional tools Several additional tools are included in the package: 

1. Dinv[j [. . .] ,sp[p,q]] — the derivative with respect to the invariant con- 
structed of the external momenta, Eq. (2.6). 

2. RaisingDRR [basis , . . . ] — the right-hand side of the dimensional recurrence 
relation j^'^"^) (basis,...) = . . . , Eq. (2.9). Note that the factor jj,^ = —1 for 
Minkovskian metrics and odd number of loops should be taken into account 
manually. 

3. LoweringDRR [basis, . . .] — the right-hand side of the dimensional recur- 
rence relation (basis, ...) = . . ., Eq. (2.10). Note that the factor jj,^ = —1 
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for Minkovskian metrics and odd number of loops should be taken into account 
manually. 

4. FeynParUF [js [basis ,...]] — the functions U and F entering the Feynman 
parametrization of the integrals in the given sector. 

Learning more There are several reduction examples in the directory Exam- 
ples of the archive file. One is encouraged to examine these examples for some hints 
of the package usage. Another good starting point to know more about the functions 
of the package is to submit a command ?LiteRed'*. 

4 Implementation notes 

The LiteRed package depends on small packages Types , Numbers, Vectors, and 
LinearFunctions, which are also included in the distributive. The Types package 
allowas one to define types and their transformation rules (e.g. vector plus vector is 
a vector). The packages Numbers, Vectors introduce specific types for number and 
vector variables. The package LinearFunctions contains the function LFDistribute 
which distributes linear functions over sum and pulls out the numbers and expressions 
having type Number. 

AnalyzeSectors uses Criterion 1 from Ref. [8] for the determination of zero sec- 
tors. FindSymmetries uses a combined approach based on the Feynman parametriza- 
tion and on the loop momenta shifts. In the first stage FindSymmetries finds map- 
pings between the simple sectors using approach based on Feynman parametrization 
very similar to the one described in Ref. [31]. This procedure works sufficiently 
fast even for complicated examples. E.g., for the four- loop onshell mass operator 
topologies the typical working time is a few minutes. 

As it was stressed above, the package performs a heuristic search of the reduction 
rules. The result of this search may strongly depend on the order in which the func- 
tions Da are listed, as well as on the choice of the irreducible numerators. Therefore, 
in case the program fails to find reduction rules, it makes sense to try changing the 
irreducible numerators, as well as the listing order of Da- 

5 Conclusion 

In this short note we have presented a Mathematica package LiteRed performing the 
IBP reduction of the multiloop integrals. The package is based on the heuristic search 
of the reduction rules (the procedure Solve j Sector) and further application of the 
rules found to the reduction problem. If the heuristic search finishes successfully, the 
rules found present a very effective solution of the reduction problem for the given 
class of the integrals. Similar to the algorithm of costruction of s-bases in Ref. [12], 
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the heuristic search of the reduction rules is not proved to terminate (and, in fact, 
seems to be not terminating in some comphcated cases). However, it appears that 
the search of the reduction rules, as implemented in LiteRed, succeeds for a larger 
class of physically interesting cases, see Examples folder in the distributive. 

The package can be useful also for some other purposes. In particular, the 
procedure FindSymmetries can be used to find mapping between equivalent sectors. 
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in[i]:= << Lit eRed ^ (it Loading the package*) 

**************** Li t eRed vl . ******************** 

Author: Roman N. Lee, Budker Institute of Nuclear Physics, Novosibirsk. 
Release Date: 12.12.2012 

LiteRed stands for Loop InTEgrals REDuction . 

The package is designed for the search and application of the 

Integration- By- Part reduction rules. It also contains some other useful tools. 
See ?LiteRed^* for a list of functions. 

in[2i:= SetDirectory[NotebookDirectory[]] ; {*Setting working directory*) 
SetDim[d] ; (*d stands for the dimensionality*) 
Declare[{l, r, p} , Vector]; {*vector variables*) 
^P[Pf P] = -^Z {* constraint*) 
in[4]:= NewBasis[p2, {sp[p-l] - 1, sp[p - 1 ~ r] - 1, sp[p ~ r] - 1, 1, r} , 

{1, r} , Directory-* "p2 dir"] ; {*Basis definition. The option Directory-* 
"p2 dir" determines the directory where all definitions for the basis will be saved*) 

Valid basis . 

Ds [p2] denominators , 

SPs[p2] scalar products involving loop momenta, 

LMs[p2] loop momenta, 

EMs[p2] external momenta, 

Toj[p2] rules to transform scalar products to denominators. 
The definitions of the basis will be saved in p2 dir 

DiskSave::dir : The directory p2 dir has been created. 

in[5]:= GenerateIBP[p2] (*IBP generation*) 

Integration-By-Part&Lorentz- Invariance identities are generated. 
IBP [p2] integration- by- part identities, 
LI[p2] Lorentz invariance identities. 

in[6]:- AnalyzeSectors[p2] {*Zero and simple sectors determination*) 

Found 15 zero sectors out of 32. 

ZeroSectors[p2] zero sectors, 
NonZeroSectors[p2] nonzero sectors, 

SimpleSectors[p2] simple sectors {no nonzero subsectors) , 

BasisSectors[p2] basis sectors {at least one immediate subsector is zero), 
ZerojRule[p2] a rule to nullify all zero j[p2 ] . 

in[7]:- FindSymmetries[p2] {*Finding unique and mapped sectors*) 

Found 8 mapped sectors and 9 unique sectors. 
UniqueSectors[p2] unique sectors . 
MappedSectors[p2] mapped sectors. 

SR[p2] [ ] symmetry relations for j[p2, ] from UniqueSectors [p2] . 
jSymmetries[p2, ] symmetry rules for the sector js[p2, ] in UniqueSectors [p2] . 
jRules[p2, ] reduction rules for j[p2, ] from MappedSectors [p2] . 

in[8]:= Timing[SolvejSector/@ UniqueSectors [p2] ] 

Sector js[p2, 0, 1, 1, 0, 0] 

Master integrals found: j[p2, 0, 1, 1, 0, 0] . 

jRules[p2, 0, 1, 1, 0, 0] reduction rules for the sector. 

MIs[p2] updated list of the masters. 



Sector js[p2, 1, 1, 1, 1, 1] 

Master integrals found: none. 

jRules[p2, 1, 1, 1, 1, 1] reduction rules for the sector. 
MIs[p2] updated list of the masters. 

Out[8]. {10.624, {1, 1, 0, 0, 1, 0, 0, 0, 0}} 

in[9]:= DiskSave[p2] ; 
QuitU 

DiskSave::overwrite : The file p2 dir/p2 has been overwritten. 

Figure 1. A simple program of finding IBP rules. 
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in[i]:= << LiteRed' (*Loading the package*) 

i^-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k Lj. fc eRed Vl . kkkkkkkkkkkkkkkkkkkk 

Author: Roman N. Lee, Budker Institute of Nuclear Physics, Novosibirsk. 
Release Date: 12.12.2012 

LiteRed stands for Loop InTEgrals REDuction. 

The package is designed for the search and application of the 

Integration- By- Part reduction rules. It also contains some other useful tools. 
See ?LiteRed'k for a list of functions. 

ln[2]:= SetDirectory[NotebookDirectory[]] ; (aSetting working directory*) 
SetDim[d] ; (*d stands for the dimensionality*) 
Declare[{l, r, p} , Vector]; {*vector variables*) 
sp[P/ p] = 1; (,* constraint*) 

ln[4]:= << "p2 dir/ p2" (*Loading basis*) 

Out[4]= p2 dir 

ln[5]:= Timing[IBPReduce[j[p2, 1, 2, 3, 4, 5], "filel"]] 

(*Reduction. Second argument is the file name to save the result*) 

Out[5]= {30.093, - [9 [-19405749861 751 770316800 + 63 337 852 789 213 794 140 160 d - 
96 838 493 039 514 249 854 976 <f + 92 253 915 806 057 780 477 952 d^ - 
61 441 859 881 785 836 896 256 d^ + 30 423 441 751 985 908 285 440 d^ - 
11 626 976 809 896282 329 088 d^ + 3 513 911 887 978 934 254 592 d^ - 
853 504 812 879 165 864 192 d^ + 168 392 515 692 831 666 432 - 
27 160529195 747 468416 d^° + 3 592160 347 775198 784 d" - 

389515 934 714 662 752 d^^ + 34 517 755 931 398 912 d" - 2 483 176 910 815 800 d'^^ + 
143 450 814 863 448 d^^ - 654 6 228 381 621 d^* + 230251 567 545 d^'' - 
6 Oil 210 070 d^^ + 109 500 174 d^^ - 1240 029 <f° + 6561 cf^) j[p2, 0, 1, 0, 1, 1]) / 
{67108 864 {-13 + d) {- 12 + d) {- 11 + d) {- 9 + d) {- 23 + 2 d) {- 21 + 2 d) {- 19 + 2 d) 
{-17+ 2d) {-15 + 2d) {-13 + 2d) {- 11 + 2 d) {- 9 + 2 d) {-7+2d)) + 
( ( - 283 807 832 792 840 965 324 800 + 948 786 748 942 911447 490 560 d - 

1 462 663235 537 746 082 583 552 <f + 1 385 351 619 026 7 60 469 495 296 - 
905 881 158 492 949 979 133 120 d^ + 435 489259 998 064 093 702 272 d^ - 
159 95 7 413 669 416 8 77 1 68 624 d* + 46033302 705 355 987 593 184 d' - 
10 555172 064 614 400 902 260 d^ + 1 949 560 973 441 170 302 024 d^ - 
291 940 360 748 911 454 721 d^° + 35 537 943 518 009 683 032 d" - 

3 513 393262 609 951175 d^^ + 280 771 095 481 687 696 d^^ - 17 972 249 382 323 578 d" + 
907 850 680 750 800 d^^ - 35 362 115 15 7 090 d^* + I 024 180145 448 d^^ - 
20 759 68 7 525 d^^ + 262 701 688 d^^ - 1561 555 d^") j[p2, 0, 1, 1, 0, 0])/ 
[536870 912 {-15 + d) {- 14 + d) {- 13 + d) ^ {- 12 + d) ^ {- 11 + d) ^ {- 10 + d) 
{-9+ d)^ {-8+ d) {-7+d) {-6+d) {-5+d) {- 4 + d) {-3+d)) + 
[81 [l 056184 729 600 - 1 829 751275 520 d+ 1 379 912 93 7 4 72 - 597241 499 904 d^ + 
164 4 75 088 128 d^ - 30 146 949360 d^ + 3 728160 432 d^ - 307 422 216 d' + 
16189 632 d^ - 492 075 d^ + 6561 d^") j[p2, 1, 1, 1, 0, 0])/ 
[268 435 456 {-14 + d) {- 12 + d) ^ {- 10 + d) {- 8 + d) {-6+d) ( - 4 + d) ) } 

ln[6]:= Quit[] 

Figure 2. A typical program of applying IBP rules. 
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